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We report the observation of tunneling anisotropic magnetoresistance effect (TAMR) in the epi- 
taxial metal-semiconductor system Fe/GaAs/Au. The observed two- fold anisotropy of the resistance 
can be switched by reversing the bias voltage, suggesting that the effect originates from the inter- 
ference of the spin-orbit coupling at the interfaces. Corresponding model calculations reproduce the 
experimental findings very well. 

PACS numbers: 73.43.Jn, 72.25.Dc, 73.43.Qt 



Tunneling magnetoresistance (TMR) devices consist of 
a tunneling barrier, typically an oxide, sandwiched be- 
tween two ferromagnetic layers of different coercive fields. 
Such systems find widespread use in sensor and memory 
application as they exhibit a large resistance difference 
for parallel and antiparallcl alignment of the ferromag- 
nets' magnetization |l|| . The TMR effect relies, within 
the simplest model |2j , on the different spin polarizations 
at the Fermi energy Ep in the fcrromagncts; it is absent 
if one ferromagnetic layer is replaced by a normal metal. 
Hence it came as a surprise that a spin- valve-like tunnel 
magnetoresistance was found in (Ga,Mn)As/alumina/Au 
sandwiches 0. The origin of the effect, labeled tunneling 
anisotropic magnetoresistance (TAMR), was associated 
with the anisotropic density of states in the ferromag- 
net (Ga,Mn)As. An enhanced anisotropic magnetoresis- 
tance (AMR) effect measured across a constriction in a 
(Ga,Mn)As film was ascribed to the TAMR effect, too 0. 
In both experiments the fourfold symmetry, expected if 
the (Ga,Mn)As hole density of states is involved, was 
broken and ascribed to strain in (Ga,Mn)As. 

Here we show that a TAMR effect can also be observed 
in sandwiches involving a conventional ferromagnet like 
iron. A stack of Fe, GaAs and Au, with iron grown 
epitaxially on the GaAs tunneling barrier, shows pro- 
nounced spin- valve-like signatures. We observe a uniax- 
ial anisotropy of the tunneling magnetoresistance. De- 
pending on the bias voltage the high resistance state is 
either observed for the magnetization M oriented in [110] 
or in [110] direction. We propose a theoretical model in 
which the Civ symmetry, resulting from the interference 
of Bychkov-Rashba and Dresselhaus spin-orbit interac- 
tions, is transferred to the tunneling probability, giving 
rise to the observed two-fold symmetry. 

A sketch of the system is shown in Fig. Ufa). The 
13 nm thick epitaxial iron layer was grown on an 8 
nm thin GaAs (001) barrier by transferring the freshly 
grown GaAs heterojunction from the molecular beam 
epitaxy chamber to a magnetron sputtering system with- 
out breaking the ultrahigh vacuum (UHV). The quality of 
the interface of a sample from the same wafer was checked 



by high-resolution transmission electron microscopy [f|. 
The Fe layer was covered by 50 nm cobalt and 100 nm 
gold which serves as back contact. The wafer then was 
glued upside down to another substrate and the original 
substrate was etched away. Finally, the circular, 150 nm 
thick top gold contact was made by employing optical 
lithography, selective etching of AlGaAs, and UHV mag- 
netron sputtering. At the Fe/GaAs and the Au/GaAs 
interfaces Schottky barriers form. The barrier heights 
can depend on the preparation technics Q and were as- 
sumed to be 0.75 eV on each side, which was found for the 
Fe/GaAs interface Q- Hence the GaAs layer constitutes 
a nearly rectangular barrier allowing, e.g., observation of 
the TMR 0, Q- In total, four batches of samples which 
differ in the preparation of the Au layer (with and with- 
out H + -plasma etching step, see e.g. Ref. Jj|) or in an 
additional annealing step (150° Celsius for 1 hour) were 
investigated. As the described features are essentially in- 
dependent of these details we focus on the results of one 
sample (annealed without H + -plasma etching) below. 

The measurements were carried out at 4.2 K using 
a variable temperature insert of a 4 Hc-cryostat with a 
superconducting coil to generate the external magnetic 
field B. We used a Semiconductor Analyzer HP 4155A to 
probe the resistance drop across the GaAs barrier in four- 
point configuration. Therefore the top Au contact was 
grounded. To vary the direction M of iron, the sample 
was mounted in a rotatablc sample holder enabling a 360° 
in-plane rotation of B. The direction of B is given by its 
angle with respect to the hard [110] direction (nomencla- 
ture with respect to GaAs crystallographic directions). 
The I-V-characteristics, measured between top gold and 
bottom Fe contact, is strongly nonlinear (not shown). 
This suggests that electron transport through the bar- 
rier is, as in previous TMR experiments, dominated by 
quantum mechanical tunneling Q ■ 

Our Letter is about the anisotropy of the tunneling re- 
sistance with respect to the in-plane magnetization M 
of the iron contact. Epitaxial iron has both a cubic 
anisotropy of bulk iron as well as an uniaxial contribution 
stemming from the interface. Magnetization reversal for 
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FIG. 1: (a) Schematic picture of a Fe/GaAs/Au tunnel struc- 
ture; (b) schematic top view on the iron layer; (c)-(h) tunnel- 
ing resistance depending on the magnetic field B swept under 
different angles <j> measured at 4.2 K and -90 mV bias (Au 
contact grounded). The double step switching mechanism is 
illustrated in (b) and (e) (see text). 



an in-plane magnetic field typically takes place in two 
steps explained by nucleation and propagation of 90° do- 
main walls @. Figs. 0Jc) - IUli) display the tunneling 
resistance as a function of magnetic field B swept in dif- 
ferent in-plane directions at a bias voltage of -90 mV and 
a temperature of T = 4.2 K. Fig.QJc) shows the resistance 
for the magnetic field swept at an angle of (f> = 90°([TlO] 
direction) from negative saturation to positive saturation 
and back. The figure focuses on the interesting region be- 
tween -0.25 T and +0.25 T. A clear spin-valve like signal 
characterized by one switching event (one jump in R) 
is observed for the resistance if B is applied along this 
hard direction. If B is applied 11° off the hard [TlO] 
axis the characteristic second switching process occurs 
at ~ 0.12 T, as is manifested in Fig. [ljd). Decreasing <f> 
the second switching point is shifted towards smaller B 
[Fig. QJe)]. This two step switching process is described 
in more detail for <f> = 68° in Fig. HJe). Starting close 
to saturation at -B [point 1 in Fig. [JJb)], the average 
magnetization direction moves towards the hard mag- 
netic [110]-axis (point 2) if B is reversed and increased 
towards positive field values. In the first step M switches 
from near the easy axis closest to the original direction 
of B beyond the easy axis located 90° sideways from this 
one (point 3). Increasing B further drives M towards 



the [TlO] direction (point 4) until, in the second switch- 
ing event. M jumps near the easy direction closest to the 
new B-direction (point 5). The signal disappears if B 
is swept along an easy direction - in the present sample 
lying at <f> = 34° - [Fig. QJf )] and changes sign for B close 
to [Fig. Dig)] or along the hard [110] direction [Fig.^h)]. 

Though reminiscent of the AMR effect the results pre- 
sented here cannot be explained by the conventional 
AMR effect of the iron layer. The resistance change 
caused by the AMR effect of the iron layer of only about 
4 mil is much smaller than the the observed change in 
the tunneling resistance of about 3.5 fi. So the AMR 
effect can be excluded as physical origin of the measure- 
ment and the question for the origin of the anisotropic 
resistance remains. 

The symmetry of the anisotropic tunneling magnctorc- 
sistance becomes more explicit at higher B where M is 
forced to follow the direction of the externally applied 
magnetic field. The data displayed in a polar plot in 
Fig. [2t[a), normalized to the resistance in [110] direction, 
were taken at B = 0.5 T at a bias voltage of - 90 mV and 
T = 4.2 K. An uniaxial anisotropy evincing the shape of 
a "horizontal 8" is clearly manifested. The resistance in 
[110] direction is typically ~ 0.4% smaller than in [110] 
direction. This anisotropy of the resistance explains the 
resistance jumps observed in Figs. [ljc)-[llh): The ac- 
tual position of the (average) magnetization determines 
the resistance. The direction highlighted by triangles in 
Fig. [2ja) correspond to the directions, taken up by the 
magnetization M in Fig. [He) for the marked B values. 
The thin red line is the result of a model calculation with 
one adjustable parameter as pointed out below. 

The anisotropy depends on the applied bias voltage. If 
the bias voltage is reversed from -90 mV to +90 mV the 
"8" is rotated by 90° as shown in Fig. [Zfb). The bias 
dependencies of the resistances' angular characteristics 
is summarized in Fig. Etc). While for bias voltages V 
> 50 mV the resistance is larger for the [110] directions, 
for V < 50 mV the resistances are largest for the [110] 
directions. Similar behavior was found for all samples 
investigated. 

We now introduce a model ascribing the observed 
anisotropy to anisotropic spin orbit interaction. It has 
already been stated in reference to GaMnAs junctions 
[ESEIGjI that the effect is due to spin-orbit coupling. 
However, to capture the Gi v symmetry of the observed 
TAMR, a uniaxial strain has been invoked Q lowering 
the four-fold symmetry of the GaMnAs density of states. 



On the other hand, ab-initio calculations in CoPt [11 1 
suggest that strain is not necessary to have anisotropic 
electronic structure in layered systems. What then leads 
to the twofold symmetry of the TAMR? We argue here 
that TAMR in epitaxial systems does not need an ad hoc 
anisotropic density of states. Instead, the tunneling prob- 
ability itself is strongly anisotropic due to the interfacial 
spin-orbit coupling. We propose that the two-fold sym- 



3 



(a) 



[110] 





1 .uuu- 




0.998- 


o 
*r 
r- 


0.996- 








0.996- 




0.998- 




1.000- 


(b) 


1 .002 -, 




1.001 - 


© 

T— 






1.000 - 








1.001 - 




1 .002 - 




Experiment 
Theory 



[110] 



Vbias = 

-90 mV 



o Experiment 
^—Theory 



[110] 



Vbias = 

+90 mV 



(c) [110] 



1.002 



[110] 



[110] 




en 

Cg 0.998 
0.996 



-90° -60° -30° 0° 30° 60° 90° 

FIG. 2: a) 0-scan of the tunneling resistance at 4.2 K and 
-90 mV bias in a saturation magnetic field |B|= 0.5 T and 
|B|= fO T and a theoretical fit; b) 0-scan at +90 mV; c) 
0-scans for different bias voltages. Symbols correspond to 
experimental results for -90 mV, -50 mV, 50 mV, 90 mV and 
135 mV bias; solid lines correspond to theoretical results with 
a t = 42.3 eVA 2 ,a, = 45.8 eVA 2 ,a, = -0.6 eVA 2 ,a, = 



-17.4 eV A 



Oil 



-25.1 eV A respectively. 



mctry of the TAMR is a consequence of the anisotropic 
spin-orbit interaction (SOI) that reflects the bulk and 
structure inversion asymmetries of our system. Indeed, 
the combination of bulk inversion asymmetry (Dressel- 
haus SOI) [H, [l3| and structure inversion asymme- 
try (Bychkov-Rashba SOI) 0] in GaAs-like semiconduc- 
tor heterostructures leads to a spin-orbit interaction with 
C2v symmetry. Based on this observation we consider 
the following model Hamiltonian for describing the tun- 



nelling across our metal/semiconductor heterojunction: 
H = H + H Z + H BR + H D . (1) 



Here 



H a = -— V 



m(z) 



+ V Z , 



(2) 



with m(z) the electron effective mass [in terms of the bare 
electron mass mo we assume m = m c = 0.067 mo in the 
central (GaAs) region and m = mi = m r ps tuq in the 
left (Fe) and right (Au) regions] and V(z) the conduc- 
tion band profile defining the potential barrier along the 
growth direction (z) of the heterostructurc [see Fig.[3]Ja)]. 

The Zceman spin splitting due to the exchange field 
(in the Fe region) and the external magnetic field in the 
Fe and Au (the Zeeman energy in GaAs is much smaller 
than all the other energy scales characterizing the system 
and we can therefore neglect its effect) is given by 



A(z) 
Hz = ■ 



(3) 



Here A(z) represents the Zeeman energy in the different 
regions, er is a vector whose components are the Pauli 
matrices, and n is a unit vector defining the spin quan- 
tization axis determined by the in-plane magnetization 
direction in Fe. 

The Bychkov-Rashba SOI due to the structure inver- 
sion asymmetry at the interfaces can be written as (la| 



HbR = r J! a A a xVy ~ (JyPx)S{z - Zi), 



(4) 



where, ai (a r ) denotes the SOI strength at the left (right) 
interface zi = (z r = d). We note that inside the GaAs 
barrier, away from the interfaces, there is also a Bychkov- 
Rashba SOI contribution induced by the applied bias. 
However, this contribution is negligible for our system 
and we neglect it. 

The Dresselhaus SOI resulting from the bulk inversion 
asymmetry in GaAs is incorporated in the model through 

sin, 



the term LL3, 14 



Hd = \{OxPx 

n 



°^z ( l{z) Yz 



(5) 



where the Dresselhaus parameter 7 sw 24 eV A 3 in the 
GaAs region 

0,(13, [l7l [l8| and 7 = elsewhere. 
The current flowing along the heterojunction is given 

by 



I = 



I dEd 2 k\\T a (E, k\\)[fi(E) - f r (E)], 

•=-1.1 

(6) 

where k|| is the in-plane wave vector and fi(E) and f r (E) 
are the electron Fcrmi-Dirac distributions with chemical 
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potentials /i; and fi r in the left and right leads, respec- 
tively. The particle transmissivity T a {E, ky) is found, 
as usually, after solving for the scattering states in the 
different regions. 

Calculations for the dependence of the resistance on 
the angle 9 between the magnetization in Fe and the 
[100] direction (note that 9 = tp+ir/A) were carried out at 
zero temperature and a barrier high (measured from the 
Fermi energy) of 0.75 eV. For the Fe layer we assume a 
Stoner model with the majority and minority spin chan- 
nels having Fermi momenta = 1-05 x 10 8 cm -1 and 
kpi = 0.44 x 10 s cm -1 respectively. The Fermi mo- 
mentum in Au was assumed to be Kf = 1.2 x 10 8 cm" 1 

The values of the Bychkov-Rashba parameters ai, a r 
[see Eq. (|4"])] are not known for metal-semiconductor in- 
terfaces. Due to the complexity of the problem, a theo- 
retical estimation of such parameters requires first prin- 
ciple calculations including the band structure details of 
the involved materials, which is beyond the scope of the 
present paper. Here we assume on and a r as phenomeno- 
logical parameters. We have found that due to the large 
exchange splitting in the left (Fe) region, the calculated 
TAMR is dominated by a; ; the dependence on a r is neg- 
ligible and we can set a r = 0. This leaves only ai as a 
fitting parameter in the comparison of the theoretical and 
experimental value of the ratio iZruoi /R[uo] ■ Such a com- 
parison is displayed in Fig. [2(a) for the case of an external 
bias Vq = — 90 meV and low magnetic field B = 0.5 T. 
The agreement between theory and experiment is indeed 
very satisfactory, considering that we fit the value of ai 
(the fit is 42.3 cV A ) only for the direction <p = n/2 
- this is enough for our theoretical model to reproduce 
the complete angular dependence of R(4>)/R[uo]. Pre- 
liminary ab-initio calculations confirm qualitatively the 
above picture [2l[ . 

We have performed the same fitting procedure for 
other values of the applied voltage. The results are shown 
in Figs. [2(b) and (c), where the good agreement between 
theory and experiment is apparent. Different values of ai 
are obtained when varying the bias voltage, suggesting 
that the interface Bychkov-Rashba parameters are volt- 
age dependent (unlike 7, which is a material parameter), 
as found in other systems The interface Bychkov- 
Rashba parameter, ai, in our system changes sign at a 
bias slightly below 50 mV. 

The robustness of the fit points to the following phe- 
nomenological model of the TAMR. Averaging the SOI 
Hsoi = Hbr + Hd [see Eqs. (j4} and ((5])] over the states 
of the system one obtains H$oi ~ w (k||) ■ cr [l7|, where 
w(ky) = {aky — jkxj—akx + yk y ,0). Here a and 7 
are effective Bychkov-Rashba and Dresselhaus parame- 
ters that measures the SOI induced spin precession of 
the electrons during tunneling. There are only two pref- 
erential directions defined by n and w(kii ) [see Fig. [3(b)]. 
Therefore, the anisotropy of a scalar quantity such as the 




FIG. 3: (a) Schematics of the conduction band profile along 
the growth direction of the heterostructure. (b) A spacial 
view of the model system. The vector n determines the mag- 
netization direction in Fe with respect to the [100] direction 
(x axis). The SOI induced spin precession of the electrons 
during tunneling is characterized by the vector w(ky) (see 
text). 



total transmissivity is obtained as a perturbative series 
of n-w(k||), since the SOI is much smaller than the other 
relevant energy scales in the system. Averaging over the 
in-plane momenta to get the full current, the anisotropy 
is determined, up to the second order, by ([n ■ w(k||)] 2 ) 
[the first-order term vanishes, since w(kii) = — w(— kn)] 
. Thus, the tunneling current anisotropy is proportional 



to S7 sin 2(9 [22j. Taking into account that 9 = 4> + ir/A 



and that the observed anisotropy is small one obtains 
for the TAMR, R(cj>)/R [lw] - 1 ~ 57(cos2</> - 1). This 
is precisely the kind of angular dependence experimen- 
tally found (see Fig. [2]). Assuming that the spin-orbit 
parameters are voltage dependent, one can change the 
sign and magnitude of the anisotropy, cry, by varying the 
bias voltage, as shown in Fig. Notably, if cry w 0, one 
obtains a suppression of the TAMR effect, a situation 
corresponding to a bias voltage of 50 mV [Fig. [2(c)]. 

In summery we observed TAMR in an epitaxial fer- 
romagnet/insulator system and propose that the effect 
occurs whenever both Rashba and Dresselhaus SOI arc 
involved in the tunneling process. 

Financial support by German Science Foundation 
(DFG) via SFB 689 and by BMBF (nanoQUIT) is grate- 
fully acknowledged. 



[1] 

[2] 
[3] 

[5] 
[6] 



Electronic address: juergen.moser@physik.uni-regensburg.de 
I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

M. Julliere, Physics Letters 54, 225 (1975). 

C. Gould et ai, Phys. Rev. Lett. 93, 117203 (2004). 

A. D. Giddings et ai, Phys. Rev. Lett. 94, 127202 (2005). 

J. Moser et ai, Appl. Phys. Lett. 89, 162106 (2006). 

Y. G. Wang and S. Ashok, J. Appl. Phys. 75, 2447 

(1994). 



5 



[7] S. Kreuzer et ai, Appl. Phys. Lett. 80, 4582 (2002). 
[8] R. P. Cowburn et ai, J. Appl. Phys. 78, 7210 (1995). 
[9] C. Riister et ai, Phys. Rev. Lett. 94, 027203 (2005). 
[10] H. Saito, S. Yuasa, and K. Ando, Phys. Rev. Lett. 95, 
086604 (2005); 

[11] A. B. Shick, F. Maca, J. Masek, and T. Jungwirth, Phys. 

Rev. B 73, 024418 (2006); 
[12] G. Dresselhaus, Phys. Rev. 100, 580 (1955). 
[13] U. Rossler and J. Kainz, Solid State Commun. 121, 313 

(2002). 

[14] R. Winkler, Springer tracts in Modern Physics, Vol. 191 

(Springer, Berlin 2003). 
[15] E. A. de Andrada e Silva, G. C. La Rocca, and F. Bassani, 



Phys. Rev. B. 50, 8523 (1994); ibid. 55, 16293 (1997). 
[16] V. I. Perel et ai, Phys. Rev. B 67, 201304(R) (2003). 
[17] S. D. Ganichev et al, Phys. Rev. Lett. 92, 256601 (2004). 
[18] L. G. Wang, W. Yang, K. Chang, and K. S. Chan, Phys. 

Rev. B 72, 153314 (2005). 
[19] J. Wang, D. Y. Xing, and H. B. Sun, J. Phys.: Condens. 

Matter 15, 4841 (2003). 
[20] N. W. Ashcroft and N. D. Mermin, Solid State Physics 

(Saunders College, Philadelphia 1988). 
[21] V. Popescu and H. Ebert (private communication). 
[22] A. Matos-Abiague and J. Fabian (unpublished). 



